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EQUATIONS  OF  A SIMPLE  FLAME  SOLVED  BY 
SUCCESSIVE  APPROXIMATIONS  TO  THE  SOLUTION 
OF  AN  INTEGRAL  EQUATION  * 

G.  Klein 


ABSTRACT 


The  problem  of  an  idealized  flame  whose  underlying 
chemical  reaction  is  unimolecnlar,  reversible,  and  of  the 
first  order,  which  has  already  been  treated  and  solved  in  the 
references  quoted  below, is  reconsidered  here  (kinetic  energy 
of  the  gas  stream  being  neglected).  Its  solution  is  made  to 
depend  on  the  solution  of  an  integral  equation  which  contains 
an  unknown  parameter  whose  eigenvalue  has  to  be  determined. 
This  equation  is  solved  by  a method  of  successive  approxima- 
tions . 


* This  work  was  carried  out  at  the  University  of  Wisconsin  Naval 
Research  Laboratory  under  Contract  N7  onr-28511  with  the  Office 
of  Naval  Research. 
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1.  INTRODUCTION 

The  problem  of  an  idealized  flame  -whose  underlying 
chemical  reaction  is  uniruolecular,  reversible,  and  of  the  first 
order,  which  has  already  been  treated  and  solved  in  the  refer- 
ences quoted  below,  is  reconsidered  here  (kinetic  energy  of  the 
gas  stream  being  neglected).  Its  solution  is  made  to  depend  on 
the  solution  of  an  integral  equation  which  contains  an  unknown 
parameter  whose  eigenvalue  has  to  be  determined.  This  equa- 
tion is  solved  by  c.  method  of  successive  approximations . 

Except  for  minor  and  obvious  deviations*,  the  nota- 
tion is  the  same  as  that  used  in  the  first  reference  quoted;  equa- 
tions there  are  referred  to  on  the  left  margin. 


♦The  only  ones  being  (cf.  equations  1.  1-1.3) 


d d 

(11.  7-27)  ffr-  £ 


(-eliminating  the  distance  variables) 


(11.7-25)  * (-an  essentially  positive  quantity) 


¥ I 

(11.7-31)  b = ” ~ (-for  conciseness) 
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Flame  equations.  These  are,  in  terms  of  dimensionless  variables  and 
parameters,  the  equations  of  continuity  ( or  chemical  reaction),  diffusion, 
and  energy  (or  thermal  conduction);. 

d Q i 

(11. 1-25)  9 dr  * ~ p-  ' 1.1 

(11.7-28)  9 JF  = f (X~G)  ' 1-2 

(11. 7-31  f s ■£(&-<*„)- 1.3 


Hot  Boundary  conditions.  At  the  hot  boundary  chemical  reaction, 
diffusion  and  thermal  conduction  cease.  Thus 

(11.7-32)  R ( ,*•.)■*  O , 1.4 

(11.7-34)  Gao  15 

Equation  1.3  is  the  integrated  energy  balance  equation,  and  by  suitable 
choice  of  the  constant  of  integration  the  third  boundary  condition,  that 
the  temperature  gradient  must  vanish, 

f(c*o)-0  , 1.6 

has  already  been  taken  care  of. 

Cold  boundary  conditions.  If  one  assumes  a conventional  functional 
form  for  the  reaction  rate,  where  the  latter  does  not  actually  vanish  at 
the  cold  boundary  temperature,  some  care  is  needed  in  the  stipulation  of 
the  cold  boundary  conditions . Experimentally,  however,  and  in  computa- 
tion where  in  any  case  one  confines  oneself  to  a limited  number  of  decimal 
places,  the  reaction  rate  can  be  taken  as  zero  at  and  near  the  cold  boundary 
temperature.  Thus  in  practice  there  is  no  doubt  what  the  conditions  should 
be,  they  are  analogous  to  those  at  the  hot  boundary,  viz. 

R(xo,rc)-o 

(11.7-35)  X0-  G0al  , 

f(Cc)=° 


1.7 

1.8 
1.9 


Auxiliary  quantities.  It  is  convenient  to  define  the  known  linear 
function 

x*- 


1.  10 


I 


9.  3 


b/ 


and  the  parameter 


1.  11 
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Elimination  of  the  mass  rate  of  flow.  We  consider  the  temperature 
gradient  and  the  concentration  as  the  primary  dependent  variables. 

From  1.3,  1.5,  and  1.10, 

G - x*  + b<^  y l . 12 

so  that  if  the  temperature  gradient  is  known,  the  fractional  mass  rate  of 
flaw,  G,  can  be  readily  found. 


Fundamental  simultaneous  equations.  With  1.12,  1 10  equations 

1.1,  1.2  may  be  written  in  the  form 


1.  13 


1.  14 


These  equations  have  to  be  satisfied  simultaneously,  the  solutions  being 
subject  to  the  boundary  conditions.  It  should  be  noted  that  this  is  an 
eigenvalue  problem;  the  parameter  in  1. 13  is  not  known  and  depends  on 
the  boundary  conditions. 


Special  cases . In  the  following  two  cases  the  problem  simplifies 
considerably; 

When  I , it  is  clear  from  1.  10  that  1.  14  is  satisfied  by 

* 

X.  -t  X 
b*  i 

aud  hence  the  problem  reduces  to  the  solution  of  the  single  differential 
equation 


1.  15 


1. 16 


When$»0,  equation  1.  14  gives 
(11.  7-44) 


x 

X - x -r 
b«0 


1.  17 


which  when  substituted  into  1.13  again  leads  to  a single  differential 
equation, 


%■< 


ho\'-  A, 


o«0 


1.  18 


This  latter  equation  simplifies  further  if  the  reaction  rate  is  linear  in  the 
fuel  gas  concentration. 
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Reaction  rate.  The  reaction  rate  may  be  expanded: 

SR(x,r)' 

R(x}r)  - PCx^jO+Cx  - x<y?) 


5x 


1.19 


Jx« 


It  will  now  be  assumed  that  it  is  linear  in  the  concentration.  Thus 


R(x,r)  « (x-x*,)  f(r)  + RCx^r)  / 1.20 

where  y~(r)  is  a rapidly  increasing  function  of  C(cf.  graph  1<)  and  the  first 
term  on  the  right  is  the  predominant  one.  Clearly,  also 

f?(x,r)  ' (x-x*)f(r)  + R(x*r)  . 1.21 

It  is  useful  to  define  the  function  (cf.  graph  2) 

* \ 'j-(t)  dr  . 1.22 

r 

For  numerical  illustration  we  shall  take 

^ , x ~'/c 

(11.7-29).  R C x,r)  * xe  -<i-x)e  , 1.23 


(cf.  graph  3)  which  has  already  been  used  in  the  original  treatment  of  this 
problem;  (cf.  references  1,  2)*  so  that  in  this  case  (cf.  graph  1) 


\jr(z) 


r -<\*/9  )/r 

e + e r 


1.24 


The  present  method  is,  however,  not  restricted  to  this  particular 
functional  form. 


* The  only  constants  entering  are 

r *.20  . i>*  S.SBl  !(. 2 


02 


/S=  LOW  763 


, a i 

cf.  footnote  22,  p.  770  of  reference  1 , where  D = y_(  . ( 

The  independent  variable  r is  related  to  the  temoerature  by  c * ~~Zv* 

* 15000 


2700 


5000 


I 


Graph  1 

cf . 1.20,  1.24 

and  table  1 

i 


Graph  2 
cf.  1.22 
and  tab  le  1 . 


3000 


I 


2.  METHOD  OF  SUCCESSIVE  APPROXIMATIONS 
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Outline.  The  method  of  solution  adopted  here  is  one  of  successive  approx- 
imation to.  the  temperature  gradient  and  it  is  essentially  contained  in  the  relations 

t \>-  R(x 


(V) 


. x%(i 


b +•  S 


dx<r;  \ 

d7~  h 


2.  1 


2.  2 


(cf.  1.13,  1.14),  together  with  the  boundary  conditions. 


Discussion.  The  differential  equation  1.13  may  be  turned  into  an  integral 


equation,  so  tuat  with  the  boundary  conditions  1.6,  1.  9 one  has 

j ~J—  , 


Thus 


r-% 


ir“  R(v, O 

"T 


dr  - Cr^-c) 


t 


- <^-r)  } 4-,- f — 


J 


i i 


1 


f 

R(x,r) 


dr 


z ~c 

ca? 


dr 


2.  3 


2.4 


2.  5 


2.6 


r. 


Also,  by  1.  14,  1.  10  , 

d ( x-  X*) 


h 


U-x  ) «* 


l-S 


and  from  1.  25 


R(x,r)  - R Cx*rW*'x 


2.  7 


2.8 


ID) 

Suppose  an  approximation  cj.  to  ^ is  Known.  It  is  assumed  that 
2.7  gives  a corresponding  approximation  to  x-x  , and  ihat  the  two  together 
can  then  be  used  to  give  a better  approximation, from  2.  3,  with  the  use 
of  2.8.  Successive  approximations  to  the  parameter  are  similarly  found 
from  2.6.  The  method  is  clearly  justified  if  one  obtains  increasingly  convergent 
results  as  is  the  case  in  the  numerical  example  considered.  Since  the  convert 
gence  is  fairly  rapid,  almost  any  suitable  lowest  approximation  is  found  to 
yield  unique  results;  the  choice  of  lowest  approximations  for  particular  values 
of  <5  will  be  considered  in  the  following  sections. 
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Scheme  of  successive  approximations.  In  accordance  with  this  scheme 
we  thus  have,  when  is  given, 

d .<•*>)  x i j _ I-  S 
-(x  )-  j W(X  -x  )=  - — b y 

2.9 


2.10 


RCx  t c)*  R (x*z  ) +■  ( x(*(.  X*)  r)  f 
(V+i ) r<~  ~ 


2.  11 
2.  12 


r fvj 

(v*)  (v+o  /*  *l?fx  ,r) 

t **1  dr  ~ ^«,-0 

r " 

The  ratio  in  the  integrals  of  2.  12,  2.  13  for  C~CW  is  found  by 
application  of  L'Hopital's  rule.  (cf.  graph  4) 


G raph  4 
c£.  table  3 


Graph  5 
c f . tab le  3 


ioo 
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Diagram  illustrating  the 
rapidity  of  convergence  of 
successive  approximations 
to  the  parameter  ^ 
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i.  CASE  OF  DIFFUSION  NEGLIGIBLE 


10 


(V+i) 


* ' 'i^s> 


dz 


3.  1 


Method.  In  the  case  of  3*  0 the  problem  depends  essentially  on  the 
solution  of  equation  1.18  which  with  1.  21  becomes 

t 3% )•  + • 

Hence,  cf.  2.  12,  2.  13,  1.  22,  2.  2 , 

pr-,ff(x*r)  . ) 


3.  2 


3.  3 


cv+i)  * (y+i) 

^ • V + b 


+ b 


3.4 


R(x*  O 

1 

r d f?(x\ 

L dr 

O] 

?v  J 

' " a">\ 

*-  L J 

r * 

j by  (Zm)  + 

• r ’ ! 

) 

r -1 

\ 

3.5 


Lowest  approximation.  The  following  choice  of  a lowest  ('input') 
approximation  has  been  found  to  lead  to  convergent  results;  we  take  as 
lowest  approximation  the  solution  of  3.1  with  the  last  term  neglected, 
viz.  (cf.  2.  2) 

(c,vm)  - Cr^-r)  , 


where 


That  is, 


o 

o 


(o) 


( £«,-0 


b^(re,cj 

f *<*,*..> 


c 


Coo  - Ca 


Also,  by  3.4,  3.8,  1.10, 

c«>  . . 

x 'x"+h(r~-vi^rTj  J 

which  may  also  be  written 

*<0>~  x~  = t(r,c-J  . 

■Jr  fr  r \ 

X 1 } * oc  / 


3.6 


3.  7 


3.8 


3.9 


3.10 


x - x 

O 00 
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Both  the  functions  3.8,  3.  10  have  the  character  one  'would  expect  of  the 

actual  solutions. 


For  the  calculation  of  the  following  approximation^  , one  needs 
rRC*\ol  f dR<x*r)]  /[(£•* ) 

I r V L dr  /{  ret ’{  * 


3.  11 


4.  CASE  OF  LINEAR  DECREASE  OF  FUEL  GAS 

Method.  In  the  case  of  5*1  the  outstanding  feature  iB 


X - x 


4.1 


cf.  1.15,  and  it  remains  to  solve  the  equation,  cf.  1.  16 


4.2 


Hence  the  scheme  of  successive  approximations,  2.  12,  2.  13,  simplifies 
to  z -z 

(»M)  =2 —2. 

1 ) — nr—  dr 

J f 


r* 


<*>,>  f r~  ,r  r. 

t * ^ J „ O')  dr  - z) 

r 


f 


4.  3 


4.4 


with 


'Rfx'oj 

<1*0 

» J. 


dR(»*,0| 

: L. 


cv) 


(V-0 


- / 


4.5 


Lowest  approximation.  Here  we  have  arbitrarily  taken  the  parabolic 
approximation 

• 4.6 

*'»  o 


In  view  of  2.  5 it  is  clear  that  the  constant  coefficient  on  the  right  of 


4.6  is  immaterial;  it  may  at  worst  give  an  unrealistic  value  for  * Cl> 


Clearly 


* . - 

R ( x ,r) 

r f* 

s 

* J 

r* 

dfUx  ,o 
d r 


4.7 


0 1000  1000  3000 
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5 GENERAL  CASE 

Method.  In  order  to  simulate  conditions  of  a real  flame  one  should  take 
0<  & < I , and  the  relations  2.9-13  have  to  be  used  in  their  full 
generality. 


The  solution  of  2.  9 with  2.  10  may  be  written 

,E  ' • r-A  fr  -tv,  ar 

X .X  . jbee^r  | 3 


<Vj  X \-i_  L „ i d C 


f 


QC 


5.1 


where  is  any  value,  ze  < rm<  ro»,but  most  conveniently  taken  in  the 

neighborhood  where  one  expects  the  maximum  of  9-  to  occur.  The  suc- 
cessive approximations  are  then  obtained  fx>m  5.  1,  2.  12,  and  2,  13, 

with  the  use  of  2.  11  and  .•» 

r cv)  r R(x  ,c) 

tr  L ■ i 1‘ 


fv-o 


9 

x b 


■].- 


I 


5.2 


dr 


/ + ~ 
a 


L (it  J 


5.3 


r 

• ■ 


In,  * , 1 r ,,  CV>  * . 

JR(x  >c/  . , /,  / - * ) 
+ ^ Cr^)  - 

dr 

J L 


dZ 


s 5.4 


r(x(V;c) 


? 


dR  (x*\c) 

/ 

r j 

dc 

/ 

dc 

. 

- 

5.  5 


Lowest  approximations.  One  could  again  take  4.6  or  3.6  as 
lowest  approximations. 


Since  in  the  present  treatment  we  have  already  calculated 


l^l 


, the  following  form  suggests  itself 


and 


5.6 


which  has  in  fact  been  found  to  lead  to  convergent  successive  approximations. 
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Starting  from  5.6  use  has  to  be  made  of  the  following  formulae  which 
we  state  for  sake  of  completeness; 


;ate  for  sake  ot  completeness; 


5.7 


f 4a  1 

I'tH. 


(6) 


HI 

A t 


j. 

X 


= b 


C (tJ~ 1 ] * ^(hSx0^(c- >-0 * Hs-o  t 

. Vi 

<Jr 


+ 0-i) 


z_ 


dch- 

1 J»o 

4c 


5.8 

5.9 


and  relations 


5.  3-5  for  V = 0 , 


Wtettdi 
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6.  ON  THE  NATURE  OF  THE  SOLUTION 


Equation  1.  13,  and  more  concretely,  equation  4.  2,  suggest  a general 
consideration  of  the  differential  equation  of  the  form 


where 


If  ( I- if  ')  - F 
F - F(x) 


increases  asymptotically  from  zero,  rises  to  a maximum,  and  has  its 
first  (and  possibly  only)  zero  at  Xm  x,  , so  that 


6.  1 


6.2 


(f  K,'°>  (F'J,, 

<0  . 

6.  3 

Throwing  6.  1 into  the  form 

6.  4 

one  obtains  immediately  the  following  propertie 
integral  curves; 

s connected  with  the 

Locus  of  stationary  points; 

y-  F;« 

) , ( X * *,  ) ; 

6.  5 

Locus  of  points  of  infinite  gradient; 

y-  0 

,(*  + *•>  j 

6.6 

Locus  of  points  of  gradient  unity: 

X 

« 

X 

, (y  + o)  • 

6.  7 

By  the  standard  method  one  finds  from  6.  1 
of  inflexion; 

r 

the  locus  of  the  points 

y-i  7.0 

/- 4F'  ) 

6.  8 

which  is  imaginary  where  F 1 > 


There  is  evidently  a singular  point  for  the  integral  curves  at  (x,,0)and 
it  is  found  to  be  a saddle  point.  From  6.  1 it  is  easily  shown  that  the  two 
singular  so’utions  which  we  denote  by  ^ and  ^ where 


0 


■ W,'>,  ‘0 


6.10 
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have  their  first  two  derivatives  at  ( x,  t o)  given  by 


( 


j. 

X 


( 1 ± V l - F ' 


) 


6.  11 


6.  12 


It  is  the  singular  solution  u which  is  of  interest  in  the  flame  problem. 
Cf.  Graph  10.  ® ' 


at 


Also,  the  locus  of  the  points  of  inflexion  has  its  first  two  derivatives 
( Xno)  given  by 


(y'X  - i 


^ t-  » 

■*V  r 


V ' 
/-^y'  . 


6.  13 
6.  14 


so  that  the  two  branches  of  the  locus  of  the  points  of  inflexion,  Y,  and 
r , are  tangent  to  the  singular  solutions.  One  finds  the  following 
properties: 


At  ( *,  >°) 
{-sign) 


±_ 

r" 


>o 


M 


> 0 


6.  15 


(+  sign) 


0 


7 


6.  16 
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Results 

The  validity  of  the  method  in  this  simple  case  suggests  that  even  more  general 
flame  problems  may  be  treated  as  integral  e quation  problems  to  be  solved  by  a method 
of  successive  approximations. 

In  the  particular  case  considered  the  eigenvalue  parameter  q , related  to  the 
flame  velocity,  appears  to  depend  linearly,  cr  very  nearly  so,  on  <5  -which  is  re- 
lated to  the  diffusion  coefficient.  (If  this  dependence  is  rigorous,  it  should  be 
possible  to  justify  it  analytically.) 

Generalizations 

It  may  be  of  interest  to  investigate  the  (minor)  modifications  required  when  the 
specific  heats  of  the  components  are  unequal,  and  when  the  thermal  conductivity 
and  diffusion  coefficient  are  not  constants;  also  the  case  where  the  kinetic  energy 
of  the  gas  stream  is  not  neglected. 

It  would  be  of  more  immediate  interest,  however,  to  essay  the  method  on  a more 
general  idealized  flame  problem  in  which  the  chemical  reaction  is  of  the  second  order. 
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List  of  Symbols 

(in  the  order  in  -which  they  occur) 

(11.7-27)  , reduced  temperature  gradient 

(11.7-16)  r , reduced  temperature 

P , reaction  rate,  (decrease  of  fuel)  cf.  footnote  p.  1 
(11.7-29)  -f  s reaction  rate  (production  of  fuel)"  H M 

b > a constant,  related  to  specific  heat,  cf.  H " 

y , ratio  of  specific  heats  n n h 

j6  , reduced  energy  of  reaction,  1.23  H M M 

S , fraction  of  the  mass  rate  of  flow  of  fuel 
jJL  , reduced  mass  rate  of  flow 
X , mole  fraction  of  fuel 
b , reduced  coefficient  of  diffusion 
<?  , eigenvalue  parameter,  1.11 
~>jr  , a function  of  C , 1.20,  1.24 

*5^  , the  Latter  function  integrated,  1,  22 
(11.7-16)  T • absolute  temperature,  cf.  footnote  p,  4 


(11.  7-22) 
(11.6-3  ) 
(1.1.  7-19) 

(11. 7-23) 


Suffices  0 and  oo  refer  to  the  cold  and  hot  temperature  boundary. 
Exponents  (V)  refer  to  the  order  of  approximation. 


6 only: 

X , independent  variable 
ij.  , dependent  variable 
r , a function  of  X 

X , y , ordinates  of  various  loci  - not  solutions  of  6.  1 
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TABLE  1 . 
(Basic  Functions) 


*-  4 1 

Rfx  ,z))0  I 


.000  5 
.000  1 
. 000  0 z 
.000  0 
.000  0 
.^000  0 


.000  0 


Graph  I 
p.  5 

Eq.  1.  24 


Graph  2 
p.  5 

Eq.  1.  22 


Graph  S 
p . 15 
Eq.  1.  17 


Graph  3 

p.  6 

Eq.  1.  16 


dR(x*c)1 
. dr  L 


dr  -lra 
03S  23 s 


.027  72 
. 050  36 
.068  38 
.082  25 
.092  51 

.099  48 
. 103  72 
. 105  51  M 
. 105  35 
. 103  47 

. 100  31 
.096  05 
.091  02 
.085  38 
. 079  34 

.073  02 
.066  54 
.059  94 
.053  31 
. 046  66 

.040  00 
.033  33 
. 026  67 
.020  00 
.013  33 
.006  67 


.021  61 
.040  24 
. 056  00 
.069  00 
.079  40 

.087  31 
.092  92 
.096  41 
.097  95  M 
.097  77 

.096  09 
.093  11 
.089  07 
.084  19 
.078  66 

.072  66 
.066  38 
.059  89 
.053  29 
. 046  66 

. 040  00 
.033  33 
. 026  67 
.020  00 
.013  33 
.006  67 


.004  391 


1 

3 

4 


. 700  4 

. 775  9 
. 836  2 
. 883  0 
.919  0 
.945  4 

.964  8 
.978  2 
.987  2 
.992  9 
.996  3 

.998  3 
.999  3 
.999  7 
. 999  9 
1 .000  0 P 

1 . 000  0 
1 . 000  0 
1 . 000  0 
i . 000  0 
i . 000  0 
1 .000  0 


h*  3/f 

om 

. 004  4 

0 

.021  7 

.049  5 

.006  6.  . . 

.040  4 

.094  8 

.013  3. . . 

. 056  2 

. 140  0 

. 02 

. 069  4 

. 185  0 

.026  6.  . . 

. 079  9 

. 229  7 

. 033  3.  . . 

. 087  9 

. 274  1 

.04 

. 093  5 

. 318  1 

. 046  6.  . . 

.097  1 

. 361  8 

.053  3.  . . 

. 098  5 M 

. 405  0 

. 06 

.098  3 

.447  7 

. 066  6.  . . 

.096  6 

. 489  8 

.073  3.  . . 

. 09J  5 

. 531  0 

.08 

. 089  4 

. 571  7 

. 086  6.  . . 

.084  6 

.611  7 

.093  3.  . . 

. 078  8 

.650  8 

. 10 

.072  7 

. 688  9 

. 106  6.  . . 

.066  4 

. 726  1 

. 113  3. . . 

.059  9 

. 762  0 

. 12 

.053  3 

.796  5 

. 126  6. . . 

. 046  7 

. 829  7 

. 133  3. . . 

.040  0 

. 861  3 

. 14 

.033  3 

. 891  4 

. 146  6. . . 

.026  7 

.919  3 

. 153  3.  . . 

.020  0 

.944  9 

. 16 

.013  3 

.967  8 

. 166  6. . . 

. 006  7 

. 956  8 

. 173  3.  . . 

Graph  6 

p.  12 


Graph  6 

p.  12 


Graph  8 

p.  1$ 


Graph  8 
p.  15 
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TABLE  3 

^Parameter  Values) 


Lowest 

approx. 


5*0 

Eq,  3.8 


<1 

162.9 
124.  7 
122.  5 
121.6 


901. 3 
689.  5 

677.4 
672.6 


5*  i 

Eq.  4.6 


5*  % 
Ea.  5.6 


136.0 
458.  1 

383.0 
414.4 
398.  5 
407.  2 

402.  5 

405.0 

403.  5 


752 
2 534 
2 119 
2 292 
2 204 
2 252 
2 226 
2 240 
2 232 


1 923 
1 814 
1 862 
i 849 


Values 
pre  vioualy 
obtained  by 
numerical 
integration, 
(cf,  ref.  2) 


2 240 


1 840 


Note:  Though  the  numerical  results  are  correct  to  three  significant  figures, 

no  painstaking  accuracy  has  been  aimed  at  in  the  computations  as  the  object  has 
been  solely  to  demonstrate  the  validity  of  the  method;  thus,  all  integrations 
have  been  carried  out  by  trapesoidal  summation  with  only  27  equal  intervals. 
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NOTICE:  WES  N GOVERNMENT  OR  OTHER  DRAWINGS.  SPECIFICATIONS  OR  OTHER  DATA 
ARE  USED  FO®  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A DEFINITELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  TI1E 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION #0R  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHT'S  OR  PERMISSION  TO  MAN  UF  ACT  ITU 
USE  OR  SELL  ANY  PATENTED  INVENTION  THATMAY  IN  ANY  WAY  BE  RELATED  THERETO.. 
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